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Abstract: Motivated by recent studies of intersecting D-brane systems in critical string 
theory and phenomenological AdS/QCD models, we present a detailed analysis for the 
vector and scalar fluctuations in a non-critical holographic QCD model in the high tem- 
perature phase, i.e., the chiral symmetric phase. This model is described by Nf pairs of 
D4 and D4 probe branes in a non-critical AdS§ black hole background. Focusing on the 
hydrodynamic as well as the high frequency limit, we analytically obtain spectral functions 
for vector and scalar modes on the flavor probe. Then we extract the light quark diffusion 
constant for flavor current using three different methods and find that different methods 
give the same results. We also compute the heavy quark diffusion constant for comparison 
with the light quark case. 
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1 Introduction 

Gauge/gravity duality [1-4] has become a powerful tool in studying strongly coupled non- 
Abelian gauge theories. Since its invention, it has been expected to be helpful in un- 
derstanding some mysterious behaviors of low energy QCD. Although the gravity dual of 
the realistic QCD has not yet been established till now, the Sakai-Sugimoto model [5, 6] 
has successfully realized certain aspects of low energy phenomena of QCD: the confine- 
ment /deconfinement transition has been recognized as the Hawking-Page phase transition 
[7], and the non-Abelian chiral symmetry breaking/restoring phase transition also has a 
geometrical realization [5, 8, 9]. However, this model has the undesired KK modes be- 
cause it is based on critical string theory and has the unwanted internal space. What is 
worse is that these KK modes and hadrons have the same mass scale and there is no way 
to disentangle the two scales. One effective way to overcome this serious problem is to 
consider string theory in non-critical dimensions. Some non-critical string models of QCD 
have been constructed in [10-31]. 

In this paper we consider a model of non-critical AcISq black hole solution [11, 12, 14] 
(corresponding to the deconfined phase), in which the flavor quarks can be introduced by 
adding D4 and D4 probe branes [32, 33] into this black hole background. As the non-critical 
version of the Sakai-Sugimoto model, this model has also reproduced many properties of 
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low energy QCD. Some studies about this model have been performed in the literature, 
e.g., the thermal aspects 1 have been studied in [35] and the chiral phase transition with 
external electromagnetic fields and chemical potential have been examined in [36-38]. The 
chiral phase transition can be represented by the embedding profile of flavor branes: the 
U-shape configuration stands for the chiral breaking phase while the parallel branes and 
anti-branes states the chiral symmetric phase [35]. The hadron physics related to low 
energy QCD have been intensively studied in [12, 14, 35, 39]. 

One motivation for studying such a non-critical holographic QCD model lies in that 
this model is different from the well-studied intersecting D-brane systems in critical string 
theory and some phenomenological AdS/QCD models which lack theoretical completions. 
Special features of this model are that the dilaton is a constant, which is same as the well- 
studied D3 brane geometry, and supersymmetry is completely broken, which is same as the 
Sakai-Sugimoto model. Given that studies on hadron physics and chiral phase transition in 
this model have reproduced many similar features to its critical version, we here focus on the 
non-equilibrium state physics (not far from equilibrium state) of quark sector — the retarded 
Green's function of mesonic operator and the corresponding spectral functions. We here 
mainly probe some differences and similarities between the Sakai-Sugimoto model and its 
non-critical version. Meanwhile, one byproduct of this study is to further demonstrate 
some universal features of gauge/gravity duality as well as eliminate some non-universal 
ones which are peculiar to some models. 

Specifically, we in this work study the vector/scalar fluctuations on flavor probe, calcu- 
late their correlation functions 2 to reproduce the meson spectral functions and then extract 
quark diffusion constant. We here work in the quenched approximation, i.e., we do not take 
into account the backreation of flavor probe branes on the AdS§ black hole background. 
For simplification we will focus on the chiral symmetric phase which is also referred to as 
high temperature phase. In this phase mesons should not be confined or bounded together, 
moreover, they are no longer long-lived in contrast to mesons in the chiral broken phase. 
Oppositely, we expect them to be quasi-normal modes in the black hole backgrounds. This 
phenomenon can be described as meson melting [42, 43] in chiral symmetric phase (actu- 
ally, this idea is subtle when considering the processes of chiral phase transition because 
mesons prior to transition do not correctly correspond to quasi-normal modes after tran- 
sition [44]). In general, such a study usually needs too much numerical work and, as a 
result, applicability of results depends on the precision of numerical computation. So in 
this paper we will merely consider the hydrodynamic as well as high frequency limit, both 
of which allow for analytical calculations. 

In the hydrodynamic limit, the frequency and spatial momentum are much smaller than 
the temperature of the holographic system, and therefore the system reduces to hydrody- 
namic and we can consequently probe its hydrodynamic properties. Using gauge/gravity 
duality to study hydrodynamic properties of strongly coupled non-Abelian gauge theories 
has gained great success [45-68] (for reviews see [40, 69]). In our study we will use the 

1 Thermodynamics and probe energy loss in similar non-critical AdSs plasma have been studied in [34]. 

2 Studies about retarded Green's functions under the gauge/gravity framework have gained much atten- 
tion, for reviews see [40, 41]. 
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method proposed in [45] to get the hydrodynamic limit for meson spectral functions, whose 
form is in agreement with expectation from boundary thermal field theory, and then we 
extract the quark diffusion constant from the correlation function in the diffusive channel. 
To check correctness of such an extraction, we use another two different ways to compute 
this diffusion constant and find that three different methods give same results. This con- 
firms the universality of the membrane paradigm, which is proposed in [48] and further 
developed in [61, 62] for study of transport coefficients in strongly coupled thermal field 
theories from gauge/gravity duality approach. Because the quark introduced by probe 
D-branes is light (here, precisely speaking, it is massless), we refer to this quark diffusion 
as the light quark diffusion which is different from the heavy quark case being introduced 
later. 

We will also study the high frequency limit using the WKB approximation. In this 
regime, we analytically obtain spectral functions for different modes. Their behaviors are 
not the same: for vector modes they are scaling as w 2 which is the same as previous results 
[64], while for the scalar one it scales as w 6 which has not been observed in previous studies. 
Actually, these scaling behaviors can also be extracted from dimensional analysis for dual 
operators. 

For completeness, we then calculate the heavy quark diffusion constant in this model. 
The heavy quark can be modeled by an open string stretching from the boundary to the 
bulk. This configuration has been used in studying jet quenched parameter [70-72], drag 
force [73] and energy loss [74, 75] for heavy quark moving through quark- gluon plasma 
(for reviews, see [69, 76]). We here just take into account a static string and extract the 
diffusion constant from the correlation function of the string's embedding coordinate, which 
is proposed in [77] and studied in [77-79] for J\f = 4 super- Yang-Mills plasma. The same 
calculation for the Sakai-Sugimoto model was done in [80]. Our result is similar to that 
of [77] but different from [80]. A recent elaborate study for heavy quark diffusion in a 
non-conformal background can be found in [81]. 

The remainder of this paper is organized as follows. In section 2, we set up backgrounds 
for our study, introducing the holographic model and notation conventions, calculating 
the charge susceptibility for later convenience and working out equations of motion for 
flavor vector and scalar fluctuations. Section 3 is the main part of our study, in which we 
intensively analyze the fluctuations under the hydrodynamic limit and extract the light 
quark diffusion constant in three ways. Section 4 is concentrated on the study of spectral 
functions in the high frequency limit using WKB method. In section 5 we turn to the 
calculation of the heavy quark diffusion constant. The final section gives some discussions 
and conclusions as well as some open questions. 

2 Holographic setup and notation conventions 
2.1 Non-critical AcISq black hole background 

We take N c coincident D4-branes wrapped on a circle in six dimensions. If we take the 
large N c and near horizon limits for this brane system, they should be described by clas- 
sical solutions to the type IIA supergravity in six dimensions. To reproduce the con- 
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finement/deconfinement phase transition, the color branes have been chosen to wrap one 
circle and we will impose periodic boundary conditions for bosons and antiperiodic ones 
for fermions along the wrapped circle to explicitly break supersymmetry [7]. Then the 
gravitational background is dual to five dimensional gauge theory and when energy scale 
is much smaller than inversion of the circle radius it is effectively four dimensional. Due 
to the antiperiodic boundary conditions imposed for fermions, the left freedom in the low 
energy case is just the adjoint gluons. Thus we are left with a non-supersymmetric pure 
gauge theory with gauge group SU(iV c ) in the low energy limit. To describe quarks in 
the fundamental representation of the gauge group, one introduces Nr pairs of D4 and D4 
branes into this background following [32, 33]. As mentioned in the introduction, in this 
work we will take the assumption that Nj <C N c and then the backreation of flavor probe 
branes on the AdS§ black hole background can be ignored, i.e., we do the calculations 
under the quenched approximation. As in the Sakai-Sugimoto model, the flavor probe is 
transverse to the wrapped circle and branes configuration has been summarized in Table 1. 



Table 1. Branes configuration 





t 


x 1 


x 2 


'1 4 

X X 


u 


color D4 


X 


X 


X 


X X 




flavor D4/D4 


X 


X 


X 


X 


X 



Now we directly write down the background geometry for the deconfined phase using 
the notation conventions of [35]: 

ds 2 = (|) 2 (-/(n)^ 2 + ^dx^ + dxi) + (^ 2 ^, (2.1) 

— J cIxq A dx\ A dx2 A dx% A du A dx^ (2.2) 

where i = {1,2,3}, Fq is the 6- form field strength of 5-form field coupled with the color 
branes, Q c corresponds to the number of color branes N c and <j) is the dilaton vacuum. The 
space spanned by t and u has the topology of cigar with the minimum ut at the tip. To 
avoid such a conical singularity, t has to be periodic 

AttR 2 n AttR 2 . „ 

t~i + /3~i + - , p = - . 2.4 

OUT OUT 

Then we can take T = 1/(3 as the temperature of the black hole background geometry 
and interpret it as the temperature of the dual boundary field theory. The corresponding 
background geometry of confined phase can be found in [35]. 

We now briefly describe some physical aspects of this non-critical holographic QCD 
model. The main conclusions of this paragraph can be found in [35, 39]. As mentioned 
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in section 1, there are two kinds of phase transition for this model — one is the confine- 
ment/deconfinement phase transition and the other is the chiral phase transition, which 
is same as the critical Sakai-Sugimoto model. The former transition is realized as two 
different background geometries according to Witten [7]. Referring to the second one, it is 
geometrically realized as a thermal phase transition on the dual gravity side, marked by the 
shape of embedding profile of flavor D4/D4 branes. This can be justified by comparison of 
on-shell actions (which are related to the free energies for the two phases respectively) for 
different flavor embedding profiles. The parallel D4/D4 branes stands for chiral symmetric 
phase, while the U-shape flavor profile says that chiral symmetry is breaking. This identi- 
fication can also be confirmed by studying the meson spectrum under these two different 
embedding profiles. As in the critical Sakai-Sugimoto model, one can find that the meson 
spectrum is discrete and gapped when the embedding profile is of U-shape; while for the 
parallel flavor embedding profile, there is no discrete meson spectrum at all. Generally, 
there is no parallel flavor profile in the confined phase due to U-shape topology of the 
confined background geometry in the x^ — u plane. When one increases the temperature to 
the deconfined phase, there can be two kinds of flavor embedding profiles — parallel versus 
U-shaped. Therefore, the high temperature corresponds to the chiral symmetry restor- 
ing phase while the low temperature to the chiral broken phase. One striking feature of 
non-critical holographic QCD model is that the curvature of supergravity background is 
of order one and even cannot be suppressed by taking large N c limit. This is necessary to 
avoid an unwanted gap between low and high spin hadron spectra. However, unlike the 
case of the critical supergravity there is a problem with the possible stringy corrections 
to this non-critical supergravity. One may concern that the large stringy correction to 
the non-critical supergravity may affect every calculation on the gravity side, however the 
results in [12] seem to indicate that this non-critical supergravity approximaton to the low 
energy behavior of non-critical superstring is meaningful. Another good feature of this 
non-critical model is that its dilaton vacuum is constant and does not need to be elevated 
to M-theory setup in the ultraviolet region as in the critical Sakai-Sugimoto model. 

For later convenience of extraction of the light quark diffusion constant, we now turn 
to the calculation of the charge susceptibility in the chiral symmetric phase. To achieve 
this, we turn on the chemical potential 3 and its corresponding charge density which can 
be simulated by the time component of U(l) (in our study we will not consider the non- 
Abelian part of \J(Nf) flavor gauge group) gauge field on the flavor worldvolume. So in the 
following computation, we just assume Aq(u) ^ and set other components of the flavor 
gauge field to be zero. The total effective action for D-brane includes the Dirac-Born-Infeld 
(DBI) part and the Chern-Simons part, but we here omit the Chern-Simons term due to 
the reason given in [35] 4 : if Chern-Simons term were included in the total action, free 

3 Chiral phase transition of this model at finite chemical potential has been studied in [37]). Note that 
conventions used in our study have some differences from [37]. 

4 In actual fact, the Chern-Simons term only has some effect on the embedding profile of the flavor probe 
branes in the chiral broken phase, which will not be considered in this work. However, the qualitative 
feature of the profile in the chiral broken phase is still unchanged when adding the Chern-Simons term, i.e., 
it is still of U-shape. While in the chiral symmetric phase, we will set x'±{u) = and then contribution to 
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energies difference between chiral symmetric and chiral breaking phases is divergent. Once 
omitting the Chern-Simons term, the action for flavor probe branes is: 

Sos, = -T t N f J rfWV-d-fo^W , (2.5) 

where gi n d is the induced metric on the flavor worldvolume, T4 is the tension of flavor D4- 
brane, l s is the string length and the only nonzero component of field strength is Fq u = F u q. 
The induced metric is given by 

ds 2 ind = (|) 2 (-f(u)dt 2 + SijdjdaP + x' 4 2 du 2 ) + f^j j^, (2.6) 

where the prime means derivative with respect to the radial coordinate u and x 4 {u) is the 
flavor brane embedding profile in this AdS§ black hole background. Since we will focus 
on the chiral symmetric phase, we can set x'^ to be zero and the action can be greatly 
simplified 

Sdbi = VaSdbi, (2.7) 
Sdbi = tdu(l) Z Jl- (2nP s A'o? • (2-8) 



e<? J \R 

Here V4 is the volume of the four dimensional Minkowski spacetime where the boundary 
field theory lives, and Sdbi is the action density. Now we can get the equation of motion 
for Ao(u): 

^©' 7 £^-, (2.9) 

T^±2irl 2 s d = d, (2.10) 

where we have chosen the integration constant as the physical charge density d which is 
clear from holographic dictionary [82], and the rescaled density d is for convenience. 

From above equations we obtain the relation between the chemical potential fi and the 
charge density as 

2irl 2 fi = 2irl 2 s I A' du= j du / ==. (2.11) 



From eqs. (2.10) and (2.11), we know that the physical charge density d is proportional to 
the rescaled one d and [i = is equivalent to d = 0. Therefore we can calculate the charge 
susceptibility S as 

'^V (2.12) 



dd 




dd 






At=0 


dji 


d=0 ( 



the total D-brane action from the Chern-Simons term is exactly zero. What's more, the fluctuation around 
x'^{u) = just induces a linear term in the scalar action, and therefore do not contribute to the equation 
of motion. Thus, we can conclude that omitting the Chern-Simons term in this work will not change our 
main results. 
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The two parts on the right hand side of eq. (2.12) can be easily worked out from eqs. (2.11) 
and (2.10). Now we write down the final result for the charge susceptibility: 

„_ 2T 4 iV / (2< 2 ) 2 4 

We now change the notation conventions through the radial coordinate transformation 
r = ut/u. Then the background geometry changes to the following form 

ds2 = (^) 2 ^(-^ 2 + M^ + ^l) + (7) 2 ^ ) (2-14) 

e = jV2_ i? 2 = ^, /( r ) = l-r 5 , (2.15) 
and the induced metric on the flavor worldvolume in the chiral symmetric phase is 

ds 2 md = (^) 2 i (-/(r)dt 2 + SydxW) + (7) 2 (2-16) 

Under this coordinate transformation, the horizon u = ut changes to r = 1 and the 
boundary u = 00 to r = 0. This transformation will facilitate later analytical studies. 

2.2 Equations of motion for vector and scalar fluctuations 

We now extract equations of motion for the vector and scalar fluctuations which will be 
the key ingredients for later calculations. Here we just focus on the vector U(l) part of 
the U(iVj)x XJ(Nf) gauge group for the Nf flavor brane-antibrane pairs. Then the DBI 
action can be expanded to quadratic order in the field fluctuations. It is well-known that 
the vector and scalar modes do not couple to each other when the vector background is 
not turned on just as in this paper, and then we can separately analyze their fluctuations. 

First, we consider the U(l) gauge field fluctuation and expand the DBI action to the 
quadratic order in the gauge field fluctuation: 

Sgauge = ~T 4 N f J d^xdre^yj- det (g ind + 2ttP s F) 

= ~ (2 f / d'xdrV- det g ind F a ^, (2.17) 

where a and /3 denote t, x, y, z and r. We will choose the radial gauge, i.e., A r = 0. Indices 
a and j3 are contracted using the induced metric gi n d- The equation of motion for the 
gauge filed fluctuation is as follows, 

d a (V- det g ind F a P) = 0. (2.18) 



Inserting eq. (2.18) into eq. (2.17) results in the on-shell action for the gauge field part 
which is crucial in calculating the retarded Green's functions: 

.W = _ ( 2 < 2 ) T ^f J d A xdrd a [y/- det g ind A p F a rt 

= - ^J^ ^ j d A xdr{d r [V- det g ind A,F 

+yj- det gindd^ [A U F^] } , (2.19) 
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where [i and v go from t to z. 

Now we perform Fourier transformation for the gauge potential A^. Without loss of 
generality, we will choose the spatial momentum to be along the x direction 

-[L-e^A.ik^r), (2.20) 
fc M = (-w,g > 0,0). (2.21) 

Inserting eqs. (2.20) and (2.21) into eq. (2.19), we get the final result for the on-shell action 
of the gauge field fluctuation as 

(2tt1 2 s ) 2 T A N f u 2 T f d 4 k 1 
Ogauge - 2 e <P R3 J ( 2 7r) 4 r k,r)d r Ai(k,r) 

-A t (-k,r)d r Mk,r)]\Zl, (2.22) 
where % goes from x to z. Choosing j3 = r, x, y in eq. (2.18), we obtain respectively 

uA't + qf(r)A' x = 0, (2.23) 



d r V-detg ind g rr g xx A' x - ^-detg ind g tt g xx u (ujA x + qA t ) = 0, (2.24) 



d r (V-detg md g rr g yyA' y ) - y/ - det g^g™ [y?g lt + q 2 g yy ) A y = 0. (2.25) 

The A z component obeys similar equation as A y . We clearly see that equations of motion 
for At and A x couple together. Following the approach of [50], we now define El = 
uA x + qAt and Et = E y z = ojA VtZ as the longitudinal and transversal electric fields 
respectively, then they obey the following decoupled equations of motion: 

E »_(^M + iW + 5!z|/M £i . . (2,7) 
\f{r){q 2 f{r) - uj 2 ) rj f 2 {r) 

Now we can also write the on-shell action of the gauge field fluctuation in terms of the 
longitudinal and transversal electric fields El, Et, 



{27rl 2 s ) 2 T A N f f d*k 1 [ f(r) 



'-'gauge — „ jl 



(2vr) 4 r 



-^ELd r EL 



¥f(r) -oj 2 ' 

r=l 

. (2.28) 



f(r) 

~ 2 {Eyd r E y + E z d r E z 



uj- 



r=0 



In above equations, we have rescaled the frequency and spatial momentum to make them 
dimensionless 

R 2 oj R 2 q , 

u = — u = n a ^ ' 9 = — q = n a ^ ' ( 2 - 2 9 

liT 0.87Ti Ut 0.87Ti 

where T is the temperature defined below eq. (2.4). In the following sections, we will 
mainly use these dimensionless quantities to carry out our calculations. 
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Following the same procedure, we can get similar results for the scalar mode fluctu- 
ation, i.e., the fluctuation around the parallel embedding profile x' 4 = 0. For conciseness, 
we here merely list the main equations and omit the cumbersome details. The fluctuation 
around x' A = can be described as 

ds2 = © 2 h (- f{r)dt2 + + {f)w) + {^y^d^x^ 

= ds 2 {g ind ) + ds 2 (g 2 ). (2.30) 
Expanding the action to quadratic order in the scalar fluctuation <p gives 

Scalar = -T A N f J d 4 xdre~^y / - det (g ind + g 2 ) 

= (^) 2 J d^xdr^J- det g ind ^g afi d a <pdp<p. (2.31) 

We can directly write down the equation of motion for the scalar mode 

d a (V-fetgind9 aP dp<p) = 0. (2.32) 



After Fourier transformation of the scalar mode (p 

¥>(w) = J -0^e ik ^(k^r), (2.33) 
fc M = (-w, g,0,0), (2.34) 

the on-shell action of the scalar part and the equation of motion for (p take the following 
forms: 



Scalar = / (^4 ^T<p{-k,r)dMk,r 



(2.35) 



r=0 



v + m-r)^^)" =0 - (2 - 36) 

Before closing this section, we give a brief remark about the equations of motion and 
the on-shell actions. From eqs. (2.26) and (2.36) we can find that the transversal electric 
field El and the scalar mode <p have similar behaviors, and the only difference between 
them is in 1/r versus 5/r term in their equations which will affect their solutions near the 
boundary r = 0. When we choose the spatial momentum to be zero, the equations obeyed 
by the longitudinal and transversal electric fields El,Et coincide. Now we turn to the 
on-shell actions given in eqs. (2.28) and (2.35). Once we have solved these equations of 
motion, we can insert these solutions into eqs. (2.28) and (2.35) and then we do second 
functional derivatives of the actions with respect to the boundary values of El, Et and 
ip. The results are just the retarded Green's functions if we merely retain the boundary 
value of the action. This is proposed in [45] using AdS/CFT technique to calculate the 
operator correlation function of the dual boundary thermal field theory. As pointed out in 
[45] , we should only retain the boundary value of the action in order to satisfy some basic 
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properties of the retarded Green's function. This does not mean that we just need to focus 
on the boundary behavior of these equations. The last key point under this proposition is 
that all the solutions to these equations of motion should satisfy incoming wave boundary 
conditions at the horizon, manifesting the unreflecting feature of black hole. This directly 
determines some important properties of the final results which will be explicit in the 
following sections. 

3 Hydrodynamic limit and light quark diffusion constant 
3.1 Retarded Green's function from gauge/gravity duality 

In this subsection, we first give a short summary of the basic prescription for calculating the 
Minkowskian retarded Green's function under gauge/gravity duality approach, and then 
we take the scalar operator as an example to show how to put this prescription into pratice. 
The main content of this subsection can be found in the original papers such as [40, 45], 
and we take it here just for self-consistency of our work. One important idea under the 
framework of gauge/gravity duality is the Operator/Field correspondence. Explicitly, this 
means that the operator on the dual field theory side is in one-one correspondence with 
the field in curved background geometry. Another striking feature of gauge/gravity duality 
is that this correspondence is a weak/strong duality and therefore we can turn strongly- 
coupled field theoretical problems into weakly-interacting gravitational ones. What's more, 
one can use the semi-classical approximation to deal with quantum field theory in curved 
spacetime under weak coupling limit of gravity, i.e., taking the supergravity geometry as a 
background. In conclusion, this duality tells us that the main physics for a specific system 
are now encoded in the classical equations of motion for the fields in curved spacetime. 

The prescription of [45] for computing the Minkowskian retarded Green's functions 
can be simply organized as the following two steps: 

(1) Find a solution to the equation of motion for the mode fluctuation such as eq. (3.5) 
below with the following two properties: the first one is that the solution is set to be 
1 at the boundary (i.e., one should impose Dirichlet condition for the fluctuation at the 
conformal boundary); the second one is that for timelike momenta, the solution should 
take an asymptotic expression corresponding to the incoming wave at the horizon while for 
a spacelike momenta, the solution should be regular at the horizon. 

(2) The retarded Green's function is then given by a formula such as eq. (3.8) (roughly 
speaking, the second derivative of the on-shell action with respect to the boundary value 
of the corresponding mode fluctuation). In detail, this means that when one takes second 
order derivative of the on-shell action with respect to the boundary value of the mode 
fluctuation, only the contribution from the conformal boundary has to be taken. Surface 
terms from the horizon must be dropped. This part of metric affects the retarded Green's 
functions only through the boundary condition imposed for the bulk field fluctuation. 

One important feature of this prescription, which has been pointed out in [45], is that 
the imaginary part of J-(k, z) (see below) is independent of the radial coordinate. Therefore, 
the imaginary part of the retarded Green's function can be calculated at any convenient 
value of radial coordinate. The authors of [45] have also checked their prescription by 
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computing the retarded Green's functions in theories where they are known from other 
methods. 

Now we turn to a scalar operator </>(x, z) as an example to elaborate detailed calculation 
of the Minkowskian retarded Green's function under gauge/gravity duality, which can be 
found in the original paper [40]. The action to quadratic order in scalar fluctuation takes 
the following form, 

S = K f d*x I*" dz^g~ \g zz (d z <p) 2 + g^d^cp + m 2 <j) 2 ] , (3.1) 

J J Z B 

where K is a normalization constant, m is the mass of the scalar field, zb and zh are the 
locations of boundary and horizon respectively; g zz and g^ v is the gravitational metric 

ds 2 = g zz dz 2 + g^dx^dx". (3.2) 

The linearized field equation for cf>(x,z) is 

L ,3 Z U=g~g zz d z <£) + sTdM - m 2 <j> = 0. (3.3) 



V~9 

After doing partial Fourier transformation along x M for (ft(x, z) 

d 4 k 

(27T) 

the equation of motion for the scalar mode then takes the following form 



/rr k 
-^e ik - x f k (z)Mk), (3.4) 



L -,d z {V^99 zz d z f k ) ~ {gTkrK + m 2 )(f> = 0. (3.5) 



V~9 

Under above equation of motion, the on-shell action for the scalar fluctuation now 
reduces to the surface term 



S = I 7T^M-k)HKz)Mk)\ z ZT B , (3-6) 



T(k, z) = K^g zz f_ k (z)d z f k (z), (3.7) 



here we have chosen the boundary value for fk(z) to be unit for convenience. Then the 
prescription for retarded Green's function is that we should just take the boundary value 
of on-shell action. The retarded Green's function is therefore proposed to be 

G R {k) = -2F(k,z)\ z=ZB . (3.8) 

Though eq. (3.8) is just a proposal, its validity has been checked in the original papers and 
it plays an important role in many interesting problems. 

The last important point in above calculation is that the linearized equation of motion 
for (f>(k, z) should obey incoming wave boundary condition at the horizon z = zh due to 
fully-absorbing feature of black hole, which has been pointed out at the end of last section. 

However, the classical equation of motion such as eq. (3.5) cannot be solved by analyti- 
cal method in generic case. Therefore, many works have turned to numerical derivations of 
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thermal retarded Green's function under gauge/gravity duality. In some limiting cases like 
hydrodynamic or high frequency limits, we have semi-analytical methods to solve these 
complicated second order differential equations. In detail, for hydrodynamic limit, we 
can have a double series expansion of solution in terms of dimensionless frequency Cj and 
spatial momentum q; while for high frequency case, we can use WKB approximation to 
construct semi-classical wave functions for different modes. For WKB construction of the 
solution, we should carefully study these equations of motion near its regular singularities, 
the boundary and horizon, to get asymptotic solutions near these points and match these 
solutions with WKB solutions in the regular regions. This matching will finally determine 
the boundary solutions for different modes. In the remainder of this section, we will in 
detail solve equations of motion derived in section 2 in hydrodynamic limit. We will present 
the detailed construction of WKB solutions in section 4. 

3.2 Analytical result in hydrodynamic limit 

Before detailed analysis, we need connect asymptotic solutions to eqs. (2.26), (2.27) and 
(2.36) near the boundary r = with the retarded Green's functions of dual mesonic 
operators. It is simple to show that near r = exponents of these equations are (0,2), (0,2) 
and (0,6) respectively. Then their general solutions can be written as 

Er(k,r) = A T (io,q)(l + ...) + B T (u,q) (r 2 + . . .) , (3.9) 
E L (k, r) = A L (w, q) (1 + ...)+ B L (u, q) (r 2 + . . .) , (3.10) 
<p (k, r) = A v (Cj, q) (1 + ...)+ Bp (Cj, q) (r 6 + . . .) , (3.11) 

where '. . .' indicates higher powers of r. Precisely speaking, theses solutions are not the 
general ones because the difference of each pair of exponents is an integer. Inversely, the 
general solution should have logarithmic terms, which may induce contact terms in the 
retarded Green's function as pointed out in [45]. Given that these terms are not important 
in our extraction of the quark diffusion constant, hereafter we will omit them and take above 
solutions as the general ones. Inserting these boundary solutions into on-shell actions given 
in eqs. (2.28) and (2.35), we can easily obtain the following connections for different modes 
(the proportional coefficients for different modes will be explicitly given later): 

nR ,~ B (Cj,q) , . 

G n {u,q) ~ , (3.12) 

A (u, q) 

where G R (Cj, q) is the retarded Green's function for dual mesonic operator. 

Hydrodynamic limit means U) <C T and q <C T, or equivalently Cj <C 1 and 5 « 1. 
This fact allows us to follow [46, 83] to do double series expansion for Er(k,r), Ei(k,r) 
and <p(k, r) near Cj = and q = 0, which has been intensively studied in the literature. 
Near the horizon r = 1, exponents of differential eqs. (2.26), (2.27) and (2.35) are ±iu)/5. 
We should choose incoming wave boundary conditions at the horizon for different modes 
to reflect the unreflecting characteristic of black hole. This fact means that solutions to 
eqs. (2.26), (2.27) and (2.35) near the horizon should take the following forms 

Et,l (Jfe, r -)• 1) and tp (jfe, r -)• 1) ~ /(r)"^ 5 . (3.13) 
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As a result, we can make this condition more explicitly by assuming that the general 
solutions to these equations are of the following forms 



E T {k,r) 


= f(r) 


-™l b E T {r) 


E L (k,r) 


= f(r) 


-™l*E L {r) 


<p(k,r) 


= f(r) 


-<G/{ V(r). 



(3.14) 
(3.15) 
(3.16) 



Reasonably, we will also impose regularity of Et{t), El{t) and ip(r) at the horizon r = 1. 

It is not difficult to show that these new variables Et(t), Ei(r) and (p(r) obey the 
following equations of motion: 



E'+(r) + 



2iu)r 4 f'(r 
+ 



1 



+ 



f(r) f(r) r 
5iwr 8 



u 2 r 8 



+ 



E' T (r) 
f'(r) 



1\ iu>r 



f(r) f 2 (r) ' \f(r) 

2iwr 4 0J 2 f'{r) 



E'L{r) + 
3io)r 3 

+ 



1 



+ 



X2 



f(r) 



f(r) f(r)(q 2 f(r)-Cj 2 ) r 
ioj 3 r f'(r) 



E' L {r) 



0J 2 r 8 



f(r) f 2 (r) 
„. . (2iCjr A /'(! 



+ 



f 2 (r) (q2f(r)-u 2 ) 



+ 



u 2 



P(r) 



+ 



oo 2 r 8 



f(r) 
5iiur 8 



~ ) At) 



+ 



f\r) 

m 



luor 

W) 



+ 



u 2 



P{r) 

- em 



P{r) 



E T (r) = 0, (3.17) 



E L (r) = 0, (3.18) 



<p(r) = 0. (3.19) 



Before doing double series expansions for Er(r), Ei{r) and (p(r), we now briefly comment 
on these three equations. As has been remarked in section 2.2, the transversal electric field 
mode Et obeys similar equations as the scalar mode and the same thing also happens 
here. So in the following calculations, we can just focus on vector modes and their results 
can be simply generalized to the scalar case. 

For Ej-(r), its double series expansion in terms of Co and q is 



E T {r)=F {r)+CoF l {r) + 



(3.20) 



where '. . .' represents higher orders in cl> and q which are not important in our calculations. 
Putting this expansion into eq. (3.17) and counting by order of Cj give 



F[' + 



f'(r) 

m 
m 

m 



I) 



0. 



1 



, 2ir 4 , 3ir 3 
rJ f(r) f(r) 



(3.21) 
(3.22) 



From eq. (3.21) we can get Fq = C. In obtaining this final result, we have imposed 
regularity condition at the horizon for Fq. Along this procedure, it is not hard to work out 
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the solution for eq. (3.22): 

1 - r 5 



Fi(r) = 



1 - r 
5 + 1) In 



l)ln 



r 2 + -(l + VE)r + l 



" 2 - -(-l + V^r + l 



2V10 - 2V5axctan 



4r - V5 + 1 



2\/10 + 2v / 5arctan 



2V10 - 2\/5axctan 



4r + V5 + 1 
Vl0-2\/5 / 
' 5-V5 N 
. Vl0 + 2\/5. 



^10 + 2^ / 
In 3125 - V51n(5 - V^) + V51n(5 + V5) 

5 + V5 V 



-2\/10 + 2v / 5arctan 



-\/l0- 2\/5, 



(3.23) 



To get the final form of this solution, regularity together with vanishing conditions at the 
horizon for F\ have been imposed. Now we proceed by expanding the solution for Et to 
0(r 2 ) 



Er(k,r) = C 



i 

20 



2 Jl - -j=tt + 2v / 5coth~ 1 (v / 5) -51n5 



1 + -iQr 2 



(3.24) 



Then it is straightforward to write down the retarded Green's function for the transversal 
electric field mode Et in the hydrodynamic limit 



r* (n d\- 2 x i 2 <) 2 ^ N f R 2 _^t(M 

Lr TT [oj, q) — -z X — — -ro—. — . , 



oj 2 A t (w,< 



2-nl 2 ) 2 T A N f R 



OJ 



(3.25) 
(3.26) 



We now give similar analysis for the longitudinal electric field mode El, which will be 
much more complicated. Following the original proposition in [46, 83], we should perform 
double series expansion for Ei(r) as follows: 



E L {r) = G (r) + QG^r) + ^G 2 (r) + 

UJ 



(3.27) 



As in the transversal electric field mode Et case, here '. . .' represents higher order terms in 
Cj and q. We can immediately write down the equations of motion obeyed by Go(r), G\(r) 
and G2(r), but here we slightly alter our procedure, which will avoid lengthy differential 
equations. Equations obeyed by Go(r) are 



f(r) 



Gg - -G' = 0. 



Then solution to Go(r) is very simple, given by 

G (r) = C'. 



(3.28) 
(3.29) 

(3.30) 
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Equations for G\(r) and G2(r) are given by 

7'(r) 1\ „, Mr 3 



G'{ + 
G'i + 



f{r) l) G[ + W) C '-°' 
f(r) rj f(r) 



(3.31) 
(3.32) 



Fortunately, eq. (3.31) is exactly eq. (3.22) once we take Fq = C in the latter equation. 
Then solution to eq. (3.31) is just eq. (3.23) in which we should substitute C for C, i.e., 



G 1 (r)=F 1 (r)\ c =c>. 



(3.33) 



Putting the solution for G\{r) into eq. (3.32) and imposing regularity and vanishing con- 
ditions at the horizon lead to the following solution to Gi{r): 



G 2 (r) = ^(l-r 2 ). 



(3.34) 



Just like the transversal electric field mode Et case, we can now expand Ei(k, r) to 0(r 2 ) 
E L (k,r) = C 



1 + — 
20 



1 - -^=7r + 2v / 5coth- 1 (\/5) -51n5 j Q + ~ 







i 




V ~2 


Co 








V 


Co . 





~ e 



i + 



%_(£_ 

2 Co 



H — i w 
2 



00 



(3.35) 



where we follow the discussions in [50] to omit some unimportant terms of oo. The retarded 
Green's function for the longitudinal electric field mode El can be straightforwardly ob- 
tained 



Gf L (u,q) = -2x 



(27rl 2 s ) 2 T 4 N f R 2B L (Co,q) 



e<t> Co 2 - q 2 A L (Co, 



(2tt1 2 ) T 4 N f R 



00 + 



% 7,2 



(3.36) 



Note that the above retarded Green's functions are two point functions of gauge in- 
variant Et and El defined in section 2.2. It is necessary to represent them in terms of the 
original gauge field A„. This can be accomplished by the following transformation rules 



G 



n 2 r R 



Gxx — u G L L > 



iR 

G R 



r iH — r< 
^tx — ^tx 



2(~<R 

R 

tx 

r>R 



yy " zz u ^TTi 



(3.37) 
(3.38) 
(3.39) 
(3.40) 



where the quantities like G R and G R X are exactly correlation functions of the vector mesonic 
operators dual to the flavor gauge field A^, which can be recognized as propagator of the 
vector meson. 
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We now explicitly write down all nonzero retarded Green's functions in terms of dimen- 
sional variables uj and q, and we have also revealed their dependence on the temperature 
T 

Gf t {u,q) = -MT 



uj + iDq 2 
G« x {u>,q) = -NT- 



r xx 



oj + iDq 2 ' 
G*(u>, q )=G tx (u,,q) = -jVT- 



lo + iDq 2 ' 

G* (uj, q) = G zz {uj, q) = -jVTiu. (3.41) 

The constant D appeared in the longitudinal component of the vector meson correlation 
function is given by 

D = ikf < 3 - 42 > 

which is the light quark diffusion constant and will be discussed in detail later. The common 
factor j\f is 

_ 4vr (27rl 2 ) 2 T,N f R _ 3VEQ c N f 

~ 5^ " 207rg 8 l a ' ( ' 

To get the final result for this common factor, we have substituted values for D4 brane 
tension T4, dilaton vacuum and curvature radius R. 

We now give a brief remark on these retarded Green's functions for the vector mode. 
It is well-known that correlation function for global current operator can be decomposed 
into transverse and longitudinal parts in thermal relativistic quantum field theory, which 
has been explicitly revealed in our results. Moreover, in the long-time and long-wavelength 
limit, longitudinal and transversal parts of the correlation function have a universal charac- 
teristic dictated by hydrodynamics: the transversal part is nonsingular of frequency while 
the longitudinal part has a simple pole at oj = —iDq 2 , where D is the charge diffusion con- 
stant. We have seen that this feature has also explicitly appeared in our results. This once 
again confirms the applicability of gauge/gravity duality technique to studies of strongly 
coupled thermal non-Abelian gauge theories. Hereafter we will also refer to longitudinal 
electric field mode El as the diffusive channel. We will in section 3.3 state the extractions 
of the light quark diffusion constant in detail. 

Having completed calculations for the vector mode, we can generalize these results to 
the scalar mode case. The procedure and result are quite similar to that of the transversal 
electric field mode Et- We first expand ip(r) as for Er(r) 

<p(r) = flb(r) + wfli(r) + . . . (3.44) 

Then equations obeyed by Ho(r) and H±(r) can be simply worked out: 

^ + (m~r) H ' o=0 > (3 - 45) 

Hi + ( --)h[ + ^-H> + ^-H = 0. (3.46) 

V/( r ) r J f( r ) f( r ) 
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These equations are very similar to those for the transversal electric field modes Et 
and the only difference lies in 1/r versus 5/r term, but their effect on the final solution is 
not important. We then directly write down the final solutions of these equations 5 



H (r) = D' (3.47) 

where D' is a constant and 

H^r) = (^2\l 25 + 2\/5vr + lOv 7 ^ coth^ 1 VE - 100 + 5 In 5^) 

+-iD'r 5 - -iD'r 6 + 0(r 7 ). (3.48) 
5 2 

As explained in eq. (3.11), we only need to get the solution for (p(k,r) to (D(r 6 ), which 
can be obtained by inserting eqs. (3.45) and (3.48) into eq. (3.16) 



tp(k,r) = D' 



1 + ^ ( 

4 -~ 5 

5 



2 

Comparing eq. (3.49) with eq. (3.11), it is clear that 

w 



\i<Jjr G } . (3.49) 



A^(u,q) = D', B v (u,q) = — G) . (3.50) 



Then we can write down result of correlation function for scalar mesonic operator 

G w (u, q) - 2 x __^ ( .^ _ -AA T ^ , (3.51) 
and the factor M' is 



^^(f)' (3.52) 

We see that the dependence of the scalar mode correlation function on temperature has a 
surprising T 6 scaling which is different from the case of vector mode, but its dependence 
on frequency is same as the transverse vector mode. 

Before closing this subsection, we briefly list the mesonic spectral function. It is pro- 
portional to the imaginary part of retarded Green's function and defined as 

JK (w, q) = -21m G R (u, q) . (3.53) 



5 We have imposed regularity condition at the horizon for Ho(r) and regularity as well as vanishing 
conditions at the horizon for H\(r). 
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The explicit expressions for all nonzero components of mesonic spectral function are listed 
here 

2 

ft tt (a?, q) = 27VT- 



u 2 + D 2 q A ' 
9^ (u,q)=Wtx {u,q)=2ATT- 



LO 2 + L> V ' 



<R XX (u,q) = 2NT- 



OJ 3 



2 + D 2 q A 



*yy 

itpip 

3.3 Light quark diffusion constant 



Viyy (uj,q) = UK ZZ (u),q) = 2MTu) , 



<R W (w, ?) = 2AA'T 6 w . (3.54) 



Now we can extract the light quark diffusion constant based on above results. In fact, 
transport coefficients for the stress tensor or conserved charge have been widely studied in 
the literature, but those studies mainly focus on the gluon sector, e.g., shear /bulk viscosity 
and i?-current diffusion constant. The analysis in [64] has elaborately studied the flavor 
charge transport in D3/D7 brane system. In the following, we will follow [64] to use three 
different methods to extract this constant and we will find that our result is same as that 
of the D2-brane background studied in [48]. If we take a careful look at the metric for D2- 
brane and the metric used here, we find that they have the same black factor /(r) = 1 — r 5 . 
This indicates that this black factor may have a direct relationship with the diffusion 
constant. Actually, from analysis of equations of motion in previous subsection, we also 
found the crucial role of this black factor in determining final results. 

3.3.1 From the pole of diffusive channel 

This method has been briefly described below results of longitudinal mode correlation 
function. We now in detail elaborate physical picture as well as some key formulae behind 
this way. By this method, we need to study some general properties of current-current 
correlator purely from viewpoint of finite-temperature field theory. Then these properties 
should explicitly emerge from calculations on dual gravitational side. We then extract the 
diffusion constant for conserved current by matching of these properties on both side. 

We here follow [50, 64] to present the general Lorentz index structure of thermal 
retarded Green's functions of conserved current in relativistic quantum field theory. Due 
to the presence of nonzero temperature, the Poincare symmetry of relativistic quantum 
field theory reduces to translation and rotation invariance. If the conserved current is 
denoted as Jn(x), current-current correlator in coordinate space is defined as follows 

G% (x-y) = -W (x° - y°) <[J M (x), J u (y)}), (3.55) 

where (...) symbols expectation value in a translation invariant state and 9(x° — y°) is the 
step function. Then it can be Fourier transformed into momentum space, 

G% (x-y)= [ e *(-»)G« (fc). (3.56) 
J (2vr) 
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Conservation of current J^{x) implies that current-current correlator denned in eq. (3.56) 
should satisfy following Ward identity 



WG%{k) = 0. (3.57) 

In vacuum, the temperature is zero and therefore conservation of current eq. (3.57) 
demands that current-current correlator should take the following form 

G^(k 2 ) = P, u U(k 2 ), P^ = ^_^. (3.58) 

While in thermal relativistic quantum field theory, due to presence of rotation symmetry, 
it is then convenient to split the projector P^ v into transverse and longitudinal parts, 

P - P T + P L 

^o T o = 0, P£ = 0, PTj=Si3-^ (3.59) 

pL _ p _ pT 

One can easily check that k^PT^ = k^P^ u = 0, which indicates that quantities composed 
of Pj w and P£ v can naturally satisfy Ward identity. 

In thermal field theory with rotation symmetry, current-current correlation function is 
completely determined by two scalar functions 

G%{k) = P^U T (k) + P^Il L (k). (3.60) 

When II T = H L , this formula is exactly reduced to the Lorenz-invariant form eq. (3.58). 
Actually, when the spatial momenta are chosen to be zero, then the transversal and lon- 
gitudinal parts of the current-current correlator coincide and we also come to the result 
eq. (3.58). Without loss of generality we can take the spatial momentum k to be along x\ 
direction, i.e., k^ = (— lj, q, 0, 0) as in this note. Then all nonzero components of current- 
current correlator are 

G* X2 (k) = G* xa (k) = U T (uj, q ), G*{k) = ^n L (u,q), 
Gf xi {k) = Gg t (k) = ^U L (uj,q), G* xl (k) = J^U L (u,q). (3.61) 

For a system in stable thermodynamic equilibrium, the low frequency /momentum be- 
havior of scalar functions II T and II is universal and is described by effective hydrodynamic 
theory, see for example [101]. Under hydrodynamic approximation, the main conclusion is 
that H T is nonsingular as a function of oj because it does not couple to charge density fluc- 
tuations, and correlator involved conserved charge density must exhibit a hydrodynamic 
singularity whose dispersion relation should satisfy co(q) — > when q — > 0. As a result, the 
longitudinal part of current-current correlator H L (oj,q) has a simple pole at u = —iDq 2 
and this equality can be recognized as the hydrodynamic dispersion relation, where D is 
the charge diffusion constant associated with current J^(x). 

When coming to our present calculations, one can imagine that a conserved current 
Jfi{x) (here, U(l) flavor current) couples to the flavor gauge field A^(x, r), encoded in action 
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J^A^.r = 0)J-(*), Le., the boundary value of flavor gauge field M ,,r = 0) is taken 

as the source for flavor current on boundary field theory. Then form gauge/gravity duality, 
flavor current-current correlator can be equivalently calculated as the retarded Green's 
function of A^x, r) as shown in detail in previous subsection. Therefore, combining general 
properties and results for flavor gauge field two point function in previous subsection like 
eq. (3.41), we can easily find the longitudinal part of current-current correlator has a pole 
at 

u = -iDq 2 , D = — L;, (3.62) 

where D is the flavor charge diffusion constant. Therefore, diffusion constant is just a 
byproduct of two point function of vector mesonic operator. 

3.3.2 Membrane paradigm 

This way of calculating transport coefficients was first proposed in [48] and further devel- 
oped in [61, 62]. The main advantage of this approach is that explicit formulae for various 
transport coefficients can be expressed in terms of metric components. This property is 
the result of universality of membrane paradigm in studying hydrodynamic behaviors of 
strongly coupled thermal field theories through gauge/gravity duality. Here, we follow the 
conventions of [48]. There, the authors considered small fluctuations of black brane back- 
ground and formulae for different transport coefficients were obtained from holographically 
derived Fick's law. Those formulae are for gluon sector and we can extend them to probe 
brane sector, such as flavor charge transport coefficients. This can be accomplished by 
considering perturbation of flavor field living on probe brane. Suppose that flavor gauge 
field coupled with the conserved current has the action 



d p+2 x = dtd p xdr, ^ = {t,x\r) (i = l,...,p), (3.63) 



where g e ff is the effective coupling. The conserved charge diffusion constant has the 
universal formula [48, 64] 



y/-detg 



9xxgeffV~9tt9r 



r=rj< 



~ ~ 2 

, -gtt9rrg e ff ,„„,, 

dr r=F= > (3-64) 
rT sj- det g 



where '~' indicates the induced metric on the worldvolume of probe brane and tt represents 
the apparent horizon in the induced metric. With this universal formula for transport 
coefficients, the light quark diffusion in this model can be simply obtained 6 



D 



a)V(s) /(«)(£) fk 



5 _i 

/ 7 



u) 3 1.6vrT 
R) 



(3.65) 



U=Ut 



D In our calculations, we use the untransformed metric, i.e. the u coordinate. 



- 20 - 



3.3.3 Green-Kubo formula 



Besides above two methods, one use also Green-Kubo relation to compute flavor charge 
diffusion constant. As is well-known from non-equilibrium statistical physics [102, 103], 
all kinetic coefficients can be expressed, through Green-Kubo relation, as the correlation 
functions of the corresponding currents. In more detail, this relation connects the flavor 
charge diffusion constant D, charge susceptibility E and zero-momentum limit of vector 
spectral function 9\ XlXl (ui,q) as follows [102, 103] 

DE = km ^ ( ^ q = 0) . (3.66) 

The charge susceptibility S has been computed in section 2.1 and the zero limit of vector 
spectral function can be extracted from the previous subsection, eq. (3.41). Then we have 

„ 1, 9Ww,q = 0) 2T 4 N f (2irl 2 s ) 2 ul Kfrri 1 , . 

D = - lim '-= \ 1 ± Ui — -NT = . (3.67) 

Ew-rt) 2uj i? 3 1.6vrT v ; 

We see that three different ways give same results for flavor charge diffusion constant, 
which shows that gauge/gravity is a powerful tool in studying hydrodynamic aspects of 
strongly coupled thermal field theory. It can also be recognized as a check of correctness of 
our calculations. Before concluding this section, we would like to have a short comparison 
between non-critical holographic QCD model and critical Sakai-Sugimoto model on the 
main results in this section. Because related investigations in critical Sakai-Sugimoto model 
have not been completely presented in the literature, we have in previous sections mainly 
compared our results with that of D3/D7 brane system. But we have sufficient reasons to 
believe that similar results of mesonic operator correlation functions can be emergent for 
Sakai-Sugimoto model. Concentrating on flavor charge diffusion constant, it is found in 
[85] that D = • Therefore, one lesson from these studies is that light quark diffusion 
constant is proportional to T" 1 , in this sense the non-critical holographic QCD model 
brings no new features to us. But in section 5 we will find one difference between these 
two models for heavy quark diffusion constant. 



4 High frequency: the WKB method 

In this section we discuss the high frequency behavior for retarded Green's functions and 
spectral functions of mesonic operators. Some relevant studies for N = 4 super- Yang- Mills 
plasma can be found in [58]. The main idea under high frequency limit is that one can 
construct the WKB solutions to equations of motion for different modes once transformed to 
Schrodinger type, and then matching solutions at different regimes to obtain the boundary 
behaviors of wave functions ipr l,S- The direct conclusion under present study is that 
incoming wave boundary condition at the horizon can be extended to be imposed near the 
boundary r = 0. Physically, this is not difficult to understand: in high frequency limit, 
the effective Schrodinger potential is negative, i.e., there is no potential barrier to prohibit 
propagation of flavor fields in AdS§ black hole background; moreover, in the absence of 
potential barrier, there is no method to generate reflected waves near the boundary. Hence, 
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an incoming wave boundary condition at the horizon is still valid near the boundary. In the 
remainder of this section, we will demonstrate this by detailed analysis of above Schrodinger 
type equations. 

For the sake of convenience in analyzing these intricate differential equations, we first 
convert them into Schrodinger types through following transformations 



OS 2 



E L (k,r) = x \r ( + — 5 )Vl (k,r) , (4.2) 



V(k,r) = ^_^^ s (fe,r). (4.3) 

Then eqs. (2.26), (2.27) and (2.35) can be expressed in terms of ipT(k,r), ipL(k,r) and 
V'S (&> r ) as 

„ , s -3 + 4r 2 \9r a - 2r 8 + q 2 (-1 + r 5 ) + w 2 l 

0r *,r + ^ %; 2 '- *-ifr(k,r)=0, (4.4) 

4r 2 (1 - r 5 ) 

, if 75g 4 r 8 2q 2 [-15r 3 + 40r 8 - 2q 2 (l - r 5 ) + 2w 2 ] 

^ L ( ' r) + 4 \ [ g ~2 (1 _ r5) _^f + F^yfF^r^i 

-3 + 4r 2 (9r 3 - 2r 8 + w 2 ) ' 



r 2 (1 — r 5 ) 2 



Vi(fc,r)=0, (4.5) 



„, -35 + 4r 2 llbr 3 -q 2 (l-r 5 )+£b 2 ] 

M(k,r) + \ - - \ 2 hs(k,r) = 0. (4.6) 

4r z (1 — r°) 

Note that we do not take large spatial momentum limit so that we can fix q and even 
set it to be zero when taking large uj limit. Then above three intricate Schrodinger type 
equations further reduce to the following two simple ones when choosing q = 

_ 3 + 4 r 2 ( 9r 3 _ 2r 8 + 

^t,l (fc, r) + -f- — 2 hr,L (K r) = 0, (4.7) 

4r z (1 — r J ) 

„ , -35 + 4r 2 (l5r 3 + cD 2 ) 

V>S (*, r) + — - 2 hs (k, r) = 0. (4.8) 

4r z (1 — r°) 

These two equations both have only two singularities at r = and r = 1 in the range 
of [0, 1]. Away from these two points, we can directly write down the WKB approximate 
solutions to these equations. However, this kind of solutions is meaningless when going 
close to these singularities, where one should take a careful look at the asymptotic behav- 
iors of the original equations. At last, one needs to match the asymptotic solutions near 
singularities with the WKB ones. In the following discussions, we merely state detailed 
calculations for vector mode, and directly write down final results for the scalar mode for 
brevity. 

In high frequency limit, i.e., (i>l, eq. (4.7) can be further simplified when away from 
its singular points 

^T,L (k, r) + (1 ^ r5)2 ^T,L (k, r) = 0. (4.9) 
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Now we define the canonical momentum p{r) and action s{r) 

Co i' T i' r Co 

P( r ) = 1 5' S ( r ) = / P( r ) dr = / dr ~\ 5 ( 4 - 10 ) 

1 — r J JO 1 ~~ r 

and the two WKB solutions are given by 

1 



V>1 ==cos(s(r) + 0i), 



^ 2 r= sin (s(r) + <fa) , (4.11) 

where the phases 0i and <^>2 w ih be determined by matching above WKB solutions with 
the asymptotic one near the singular point r = 0. 
Near r = 0, eq. (4.7) reduces to 

i>T,L (k, r) + (~ + Co 2 ^j Vt,l (k, r) = 0. (4.12) 

Its general solution is the linear combination of the first and second kind of Bessel functions 

if,(r ~ 0) = C(l)y/rJi {Cor) + C{2)y/rY x {Cor) , (4.13) 

where the integration constants C(l) and C(2) will be determined later. 
Near r = 1, eq. (4.7) reduces to 

l , uP 

^t,l (k, r) + ^—f^ 4>t,l (k, r) = 0, (4.14) 
whose general solution is much simpler than the previous one: 

1 zui 1 i iu) 

^{r ~ 1) = (7(3) (1 - r)5~- + C(4) (1 - r)^ + ~ . (4.15) 

From eq. (4.1) and incoming wave boundary condition imposed at the horizon, we imme- 
diately obtain 

C(4) = => ip{r ~ 1) = C(3) (1 - r)5"f . (4.16) 

To match these solutions, we need asymptotic forms of p{r) and s(r) near r = and 
r = 1 respectively: 

P( r ) = , ~ 5 ~ w ( r ->• °)> 
1 — r° 

a) a; 

Pin = e ~ — ; r (r — >■ 1), 

FW l_ r 5 5(1 _ r ) ^ ; ' 

/" r Co 

s(r) = / dr = ~ wr (r — » 0), 

s{r) = / dr f ~ - [— wln(l — r) + Coin constant] (r — >• 1). (4-17) 

Jo 1 - r 5 
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The constant in the last equation is very intricate, so we do not present its expression here. 
Using asymptotic series expansion of Bessel function at large wr: 



Ji (Cor) 



TTUjr 



3vr\ 



Y\ (Cor) ~ \l —— sin [ Cor — -7- ] , (4.18) 
V ncor \ 4 J 

eq. (4.13) has the following asymptotic behavior at high frequency limit 

^(r ~ 0) ~ C(l)^J^cos (cor - ^) + C ^\[~^ s[n " t) ' ^ 4 ' 19 ^ 

Now it is straightforward to perform matching of solutions underlying different regimes. 
With r — > parts of eq. (4.17), matching of eq. (4.13) with eq. (4.11) gives 



3vr 



(4.20) 



Moreover, with r — )■ 1 parts of eq. (4.17) one can check that 

wr —^J+ism^Cor — ^ — > constant (1 — r)^~~ . (4.21) 



1 i£j 



y 7rp(r) 

This indicates that the suitable solution near r = should take the form 

ipr,L (k, r 0) = constant [\/r J\(ur) + z-v/rli (wr)] . (4.22) 
Consequently, the solution to Et,l (k, r — > 0) is 



£t,l (A;, r ->• 0) = constant^/ ^ [v'f Ji(u)r) + iy/rYi(ur)] . (4.23) 

Expanding this solution to 0(r 2 ), we obtain 

2i 

At,l [u,q= 0) = r , 

7TUJ 

£t,l (w, 9 = 0) = ^ '— '— , (4.24) 

where 7 is Euler's constant and we have omitted terms proportional to Co In Co in B (Co, q = 0) 
which we recognize as the contact terms. In other words, we only have interest in the leading 
behavior of the retarded Green's functions. Then the nonzero components of vector mesonic 
operator retarded Green's functions are given by 

G* ( U ) = G« („) = c£ M = _ ^ + (5 + 101n2-10 7W (4 26) 

The corresponding spectral functions for vector meson are 

(00) = VKy V (u) = 9K ZZ (00) = jTVw 2 . (4.26) 



4 
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Before concluding this section, we present related main results for the scalar mode. 
The appropriate solution near the boundary r = takes the form: 



ip (k, r — > 0) = constant y g- ^y/rJz(Cbr) + iy/rY^{Cjr)\ . (4.27) 

Expanding it to 0(r 6 ) gives 



Ap (u, q = 0) =- 



16i 

(-Hi + 12«7 + 6vr - 12iln2) 



£ v (£5, g = 0) = ^ ^— , (4.28) 

where we also omit the contact terms as in the vector mode case. The result for the 
retarded Green's function of the scalar meson is 

yR ,._ 6vri + (ll + 121n2-12 7 ) / 5 \° 6 



G% («) = " ( ^ ) ( 4 - 2 9) 



and the spectral function for scalar meson is 

15625 
786432tt 5 ' 

We see that the scalar mode spectral function is proportional to w 6 while it is w 2 for vector 
mode. Such a scaling behavior in the high frequency limit can also be extracted from the 
dimensions of operators dual to and ip. This strange behavior for scalar meson signifies 
some non-universal characteristics of this model or this duality. 



5 Heavy quark diffusion 

As the final part of this paper, we now turn to compute the heavy quark diffusion constant, 
which can be recognized as a comparison with light quark diffusion presented in subsec- 
tion 3.3. In gauge/gravity setup, heavy quark is introduced as the endpoint of an open 
string stretched from the boundary to the bulk in depth. Following this picture, many 
aspects of heavy quark moving through strongly coupled quark-gluon plasma such as jet 
quenching parameter [70-72], drag force [73], parton energy loss [74] have been studied 
(for recent reviews see [69, 76]). The heavy quark diffusion constant was first computed 
in [77-79] for M = 4 super- Yang-Mills plasma under the AdS/CFT technique. It has also 
been studied in [80] for critical Sakai-Sugimoto model, but we will find that the behavior 
for heavy quark diffusion is quite different between these two models. 

Before going into the detailed extraction of the heavy quark diffusion constant, we now 
briefly state how the diffusion of a relativistic heavy quark through the quark-gluon plasma 
is described by a generalized Langevin equation. We first give the physical picture from 
the boundary thermal field theory and then express the main formulas for heavy quark 
diffusion constant in terms of the quantities on the dual gravity side. 

The following statements can be found in [81] and more detailed and fundamental 
expressions can be found in [104, 105]. Due to the presence of the strongly-coupled quark- 
gluon plasma, the trajectory of a heavy quark through the plasma is believed to be described 
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by Brownian motion. The action for an external quark moving through the plasma should 
take the following form: 

S[X(t)] = S + J dTXpWir), (5.1) 

where So is the free quark action, X^(t) is the covariant coordinate of the heavy quark and 
J- represents the plasma drag force exerting on the quark. When the interesting energy 
scale is much smaller than the quark kinetic part, the microscopic degree of freedom of the 
plasma can be treated semiclassically. Then the motion of quark in the plasma is reduced 
to a classical generalized Langevin equation for the coordinate Xi(t): 



SXi(t) 

where 



/+oo 
dT0(T)C%)X,.(i - r) + C(t), (£W(0) = A ij (t - (5.2) 
-oo 



C ij (t) = Gg ym (t) ee -i([F(t),F(0)]), A«(jt) = G% m (t) ee - l -({F(t),F(0)}). (5.3) 

The first term in the first equation of eq. (5.2) is the friction part while the second one is 
for the random motion part. 

We now introduce the correlation time r c and when the time scale is much larger than 
this correlation time, random term correlation can be approximated as follows: 

A ij (t-t') « K ij 6(t-t'). (5.4) 

We also define a friction function as 

^(t)EE^(i). (5.5) 

Therefore in the large time scale (compared to the time correlation r c ) limit, eq. (5.2) 
reduces to the following local one 

^+^'^ = m (e(m j (t')) = K ij s(t-t'). ( B .e) 

In above equations, we have also defined two constant: the diffusion constant and the 
friction constant rf 3 

poo 

K v ee lim ^G% m {u) = - linj t— G%(cj), ^= dr 7 y (r). (5.7) 

In obtaining the final expression for diffusion constant, we have used the following relation 
for a system at equilibrium with a canonical ensemble at temperature T 

G% m (oo) = -coth^lmG%(oo). (5.8) 
And Gly m (oo) is the Fourier transformation of G£ yrn (t). 
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Now we need to express these results in terms of variables on the dual gravity side. 
The only thing one needs to do is to identify the appropriate bulk field that couples to 
the boundary operator J 7 , then solve the bulk equations for this field with appropriate 
boundary conditions. As mentioned at the beginning of this section, a heavy quark under 
the viewpoint of gauge/gravity duality is described by a string with endpoint attached 
to the flavor brane and then extending into the bulk. If we choose the static gauge for 
string worldsheet (i.e., a = r and r = t), then the embedding coordinate is just the spacial 
coordinate P(r, t). Under this trailing string picture for heavy quark, the appropriate bulk 
field is nothing but the quark position, i.e. the boundary value of the string embedding 
X l (r, t). Therefore, the correlation function of the drag force can be extracted from solving 
the bulk linear fluctuations around the trailing string and using the holographic method. 
In the rest of this section, we will carry out this calculation for the non-critical holographic 
QCD model. 

The string worldsheet coordinate is chosen as r = t, a = r and we just consider a 
static string, i.e., Xi and X4 are just constants, which can be proved to be compatible 
with equation of motion from the Nambu-Goto action. The induced metric on the string 
worldsheet is 

Following the procedure in [77] , we consider the fluctuation of the static string embedding. 
For simplification, we just fluctuate the X± direction X\ — > X\ + £(i,r), and then the 
induced metric including fluctuation on the worldsheet is 

ds 2 W SF = (^) 2 ^ {-f(r)dt 2 + C 2 ^ 2 + C' 2 dr 2 + 2(t>dtdr) + (| V ^- (5.10) 

where the dot represents derivative with respect to t and the prime means derivative with 
respect to r. The Nambu-Goto action up to quadratic order in fluctuation Q is 

C 2 Ut tl 



f(r) 



f(r)C 



(5.11) 



In the above expression, we have omitted an infinite constant term. Now the equation of 
motion for fluctuation is 

C(t,r) ul d H(r) ,. A , 

V((t,r) =0. (5.12) 



r 2 /(r) R 4 dr 

As in the calculation of meson spectral functions in previous sections, we perform 
Fourier transformation for £(i,r): 

C(t,r) = J ge-^CKr). (5.13) 
Under this transformation, the on-shell action and equation of motion take the forms: 

(5.14) 



_ ui, f du fir) . 

Son-shell = / (-*, V) d r C ( W , r) 



r=0 



C>,r) + 



f'(r) 21 uj 2 

J^--\t'^) + m C(u;,r). (5.15) 



-27- 



Note that the results presented here is similar to those for the vector and scalar modes, 
especially when we turn off the spatial momentum there. As we just want to extract 
out the diffusion constant following the proposition in [77], we here merely study the low 
frequency behavior of eq. (5.15), which is parallel to previous hydrodynamic limit study. 
The exponents of eq. (5.15) near the horizon are ±zo)/5 and the incoming wave boundary 
condition imposed at the horizon indicates that the physical solution should take the form 



C(w,r)~ (1-r 5 ) ^g(r), 
here, g(r) should be regular at the horizon. Then we expand g(r) in order of Q: 

a( r ) = 9o{r) + £ffi(r) + . . . 
One can get equations obeyed by go(r) and gi(r): 



flf (r) + 



f'(r) 2 

f(r) r 

f'(r) 2 

f(r) r 



9o(r) = 0, 



2ir 4 



3ir a 



9i(r) + J^9o( r ) + J^)So(r) = 0. 



The regularity condition at the horizon r = 1 solves go(r) as 

9o(r) = c 



(5.16) 
(5.17) 

(5.18) 
(5.19) 

(5.20) 



Regularity together with vanishing conditions at the horizon for gi(r) also solve eq. (5.19) 
as 



3ic 

yu ; 200 



2\/25 - wV5tt + Sv^cosh- 1 [ - 



25 In 5 



+ hcr 3 + 0(r i ). (5.21) 



In the final solution for gi(r), we have expanded it to C(r 3 ) which is sufficient for calculating 
the diffusion constant. 

The solution for £ (w,r) has the required form as in eqs. (3.9), (3.10) and (3.11): 



ICOJ 



C(u,r) = c(l + ...) + — | 



• 3 + ...). 

The retarded Green's function for the coordinate fluctuation C (w, r) is 



-2 x 



4tt/2#4 c 



(5.22) 



(5.23) 



We can first obtain the constant k, which describes the correlation of noise in Langevin 
dynamics of stochastic motion 7 , given by 



KS _ lim ^ Im0>) = ^ r 3. 
uj^O CO 25/f 



(5.24) 



7 We will not go into details of this description for heavy quark moving through quark-gluon plasma, 
which can be found in [77]. 
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Then the heavy quark diffusion constant can be extracted from the Einstein relation for 
stochastic motion: 

2T 2 l 2 q 25 1 i i 

D s — = ^isr (5 - 25) 

This result is the same as that of [77] for N = 4 SYM plasma except for the second 
factor in eq. (5.25) which is model-dependent. But our result is different from that for 
Sakai-Sugimoto model studied in [80] which has the scaling ~ T 4 in noise correlation 
k. This manifests another difference between Sakai-Sugimoto model and its non-critical 
version. However, we see that the light and heavy quark diffusion constants both depend 
on temperature as 1/T in this non-critical holographic QCD model, which indicates some 
universal features of gauge/gravity duality. 



6 Conclusion 



In this paper we investigated some properties of the chiral symmetric phase for a non-critical 
holographic dual model of QCD, which is described by Nf pairs of D4 and D4 branes in the 
near-horizon geometry of N c coincident D4 branes in six dimensions (we just focus on probe 
limit, Nf <C N c ). In such a chiral symmetric phase, the Nf flavor pairs extend parallel 
into the black hole background. Then the induced geometry on the flavor worldvolume is 
nearly the same as the background except one dimension less. Moreover, the embedding 
profile of the flavor branes is very trivial: just one pair of parallel lines in the x^ — u plane. 
In this sense, this model is very similar to Sakai-Sugimoto model and also much simpler 
than the well-investigated D3/D7 system. So our studies for flavor gauge and scalar fields 
are similar 8 to their bulk counterparts studied in [40, 45, 46, 48-51, 58, 59, 61-63]. 

We first extracted the charge susceptibility for later convenience. Fluctuations of flavor 
gauge and scalar modes were elaborately analyzed, including their equations of motion, re- 
tarded Green's functions and spectral functions in hydrodynamic as well as high frequency 
limit. Holographic spectral function and quasi-normal mode for mesonic operator in fla- 
vored large N c gauge theory or AdS/QCD model have been studied in a huge number of 
papers [42-44, 64, 65, 81, 84-100]. Due to complicatedness of differential equations such 
as (2.26), (2.27) and (2.35), in general one cannot get analytical results and studies along 
this line need much numerical work. However, the physical picture underlying this study 
is quite general and simple. In the chiral broken phase, meson spectrum is featured with 
a discrete set of stable states [14, 35, 39] and its spectral function is characterized by a 
series of 5- function peaks as pointed out in [64]. While in the chiral symmetric phase, the 
meson spectrum should be featured by a set of quasi-normal modes in the induced black 
hole metric on the flavor branes. It is regarded that this phenomenon for meson can be 
described by holographic meson melting [42, 43], but this idea is subtle [44]. Numerical 
extraction of quasi-normal modes for Sakai-Sugimoto model can be found in [85] . We leave 
detailed investigations about transition from stable mesons to quasi-normal modes in this 
non-critical model following the approach of [86] for further work. 

8 One main difference comes from the common factor denoted as M and N' in this paper. This factor 
counts the degree of freedom, it is proportional to N% for bulk modes while Q c Nf for flavor ones. 
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Therefore, we mainly studied two limits which allow us to use some analytical methods 
to get clear results. In the hydrodynamic limit, we first got analytical results for retarded 
Green's functions of different flavor modes. The retarded Green's functions for vector 
meson operators were classified according to the rotation invariance of four dimensional 
thermal field theory into transverse and longitudinal parts as illustrated in eq. (3.41): the 
transverse part of vector correlator has no pole while the longitudinal one has a simple 
pole, which is responsible for the flavor diffusion. The scalar meson correlator has similar 
behavior as that of transverse part of vector meson: scaling linearly with frequency and has 
no pole. One apparent difference between them is their dependence on the temperature T: 
for the vector mode (both the transverse and longitudinal modes) it is linearly dependent 
on T as shown in eq. (3.41) while for the scalar meson T 6 -dependent as shown in eq. (3.51). 
These characteristics are consistent with previous studies both for bulk fields and flavor 
fields. Spectral function is then just the byproduct of retarded Green's function. Although 
these mesonic correlation functions and spectral functions have the same structure as those 
of other holographic models, either from bottom-up approach or top-down intersecting D- 
brane systems, but one main difference between them is that their scaling dependence on 
the temperature and frequency, which should be recognized as model-dependent. 

We then extracted the light quark diffusion constant from the simple pole of longitu- 
dinal part of vector mode correlator. This lowest quasi- normal frequency stands for the 
dispersion relation u) = —iDq 2 where D is just the diffusion constant. We also used an- 
other two different means to compute this constant and performed the cross check. One of 
the methods is using Green-Kubo formula and the other one is using transport coefficient 
formula derived in [48, 61, 62] under black hole membrane paradigm. The same results by 
using different methods indicate some universal features of gauge/gravity duality. Another 
limit case which allows for analytical investigation is the high frequency or zero temperature 
limit. We first inverted these second order differential equations to standard Schrodinger 
types and then used WKB approximation to solve them. We found that the spectral func- 
tions are not dependent on temperature which is just the trivial result of zero temperature 
limit and their dependence on frequency can be simply extracted from dimensional analysis 
for corresponding dual operators as pointed out in [64]. 

In the final part of this paper we focused on the heavy quark diffusion constant in this 
model. This constant plays an important role in heavy ion physics that can be studied at 
RHIC and LHC experiments. We followed the prescription in [77] to concisely calculate 
this constant and found similar behavior with previous studies [77-79] for N = 4 super- 
Yang-Mills plasma but different from [80] for Sakai-Sugimoto model. This result may signal 
one main difference between Sakai-Sugimoto model and its non-critical version. 

There are some other interesting problems which are not covered in this work but 
worth of investigation in future work. The first one is about the relationship between 
stable meson state in chiral broken phase and its quasi-normal mode counterpart in chiral 
symmetric phase. In this kind of model, the flavor branes and color branes are trans- 
versely intersecting so that the bare quark mass is exactly zero. So, introducing quark 
mass and probing its effect on spectral function or meson melting in the process of chiral 
phase transition is interesting. When turning on background for U(l) gauge field to model 
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external electromagnetic fields and chemical potential, the calculations for correlation func- 
tions are more difficult due to coupling of different flavor modes. Studies along this line 
can be found in [84, 87, 88, 90, 95, 96]. Generalizing these studies to other holographic 
systems also seems interesting. We have focused on the probe limit in this paper, but 
the supergravity background in six dimension taking into the flavor backreaction has also 
been constructed in [11], so studying meson spectral functions beyond probe limit in this 
non-critical holographic QCD model also deserves explorations. 
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